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SUMMARY 


A study was conducted to determine the feasibility of a satellite 
attitude fine -control system using the interaction of the earth's magnetic 
fieid with current-carrying coils to produce torque. The approximate 
intensity of the earth's magnetic field was determined as a function of 
the satellite coordinates. Components of the magnetic field were found 
to vary essentially sinusoidally at approximately twice orbital frequency. 
Amplitude and distortion of the sinusoidal components were a function of 
satellite orbit. 


Two systems for two-axis attitude control evolved from this study, 
one using three coils and the other using two coils. The torques devel- 
oped By the two systems differ only when the component of magnetic field 
along the tracking line is zero. For this case the two-coil system 
develops no torque whereas the three-coil system develops some effective 
torque which allows partial control. The equations which describe the 
three-coil system are complex in comparison to those of the two-coil 
system and require the measurement of all three components of the magnetic 
field as compared with only one for the two-coil case. 

Intermittent three-axis torquing can also be achieved. This torquing 

S inL«frfa«. COar \ e ° r f ° r storeVZeZn 

r , , J on wlle els. Such a system has the advantage of requiring 

no fuel aboard the satellite. H ing 


For any of these magnetic torquing schemes the power required to 
produce the magnetic moment and the weight of the coil seem reasonable. 


INTRODUCTION 


For most satellites some form of attitude control is required, and 

+ l he contro1 must he fairly precise. The fine control of 
the attitude of a near-earth satellite utilizing the earth's magnetic 
field is the subject of this report. Other methods of control studied 
at Ames Research Center have used inertia wheels and gyros. Reference 1 
discusses m detail the use of inertia wheels for finf control. 
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To utilize the earth's magnetic field for control purposes, some 
knowledge of its shape is necessary. Accordingly, there will he presented 
the results of a brief study of the variation of the components of the 
earth's magnetic field at the satellite, in vehicle coordinates. It will 
then be shown how this field can be used to point the vehicle in a spec- 
ified direction. Brief mention will be made of how the field can be used 
to dump momentum from the vehicle. It will be shown that the power and 
weight required for this system are reasonable, and that the transient 
response depends, in part, on the relative orientation of the field. . 

The requirements placed on satellite attitude control systems vary 
from one satellite to another ; however, it can be stated generally that 
a tracking line in the satellite must be maintained in alinement with 
some exter nal reference line of sight to a specified accuracy. Tracking 
of the external reference line is to be maintained in the presence of 
disturbing torques and satellite orbital motion. 

As a representative example for this report an attitude control 
system utilizing the earth's magnetic field will be analyzed to deter- 
mine whether it can meet the design requirements for the fine control 
of the attitude of the Orbiting Astronomical Observatory (0A0) which 
is an earth satellite currently under study by NAfJA. This satellite 
vill carry a rigidly mounted telescope whose optical axis defines the 
tracking line. The telescope is to be used to observe stars for various 
experiments; the external reference line is therefore defined as a line 
from the satellite to the star. Precise attitude control of this satel- 
lite is required about two axes, while roll about the tracking line should 
be limited to ±^5°. The most stringent requirement of the 0A0 is a 
pointing accuracy of l/lO second of arc. 


NOMENCLATURE 


All units are in the rationalized MKS system. 

A ang ular error: the angle between the wracking line and the 

external reference line, radians 

a area enclosed by coil, meters 2 

B magnetic induction, Weber s/me ter 2 

b coil cross-sectional area, meters 2 

D coil diameter, meters 

G(s) open-loop transfer function, volts/radian 

H angular momentum, kilogram-meter 2 /sec 
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coil current, amperes 

unit vectors along the x,y, and z axes, respectively 

moment of inertia about vehicle axis, kilogram-meter 2 

sensor gain, volts/radian 

amplifier gain, volts /volt 

rate gain, volts/volt/sec 

line of sight to target 

coil circumference, meters 

magnetic moment of coil; ampere-meter 2 

mass, kilograms 

number of coil turns 

power, watts 

earth equivalent magnetic dipole moment, ampere-meter 2 
resistance, ohms 

distance from dipole center to field point, meters 

Laplace transform variable 

torque , Newton-meters 

tracking line of vehicle 

disturbance torque. Newton-meters 

orthogonal vehicle axes, x along tracking line 

dimensionless damping ratio 

permeability of free space; henrys/meter 

r e s i st i vity ; ohm-met er s 

dens ity; kilogram/ met er 3 

angle that r makes with the dipole axis, radians 
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w wheel angular velocity; radians/sec 

u> n undamped natural frequency; radians/sec 

Subscripts 


Lim 

sensor output limited 

r 

radial component 

T 

tangential component 

x, y,z 

axes of the vehicle 

o 

vector 


THE EARTH'S MAGNETIC FIELD 



The analysis of a 
cor.trol system making use 
of the earth's magnetic 
field requires knowledge of 
the field in vehicle coordi- 
nates as a function of the 
satellite's position in orbit. 
Accordingly, the equations 
of the earth's magnetic field 
at the satellite were derived 
as functions of satellite 
position. It was assumed in 
th:.s analysis that the earth's 
magnetic field can be 
expressed in terms of a 
magnetic dipole set along 
the geomagnetic axis, positive 
toward the geographic south 
pole. The declination of 
the axis was assumed to he 
17°. The satellite position 
at any instant is expressed 
in terms of the polar coordi- 
nates r and as shown in 
the adjacent sketch . 


Sketch (a) 


4r- > 
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Then the equations 
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define the components of the magnetic field at the satellite (ref. 2, 
p. 425). These assumptions give a magnetic field which is a reasonable 
approximation of the actual field at satellite altitudes (ref. 3). 

Equations (l) , plus those representing the motion of the earth and 
the satellite, must be analyzed to determine the variation of the field 
at the satellite in vehicle coordinates. Appendix A gives in detail all 
these equations, and a method for combining them to give the desired 
results. However, it is difficult to visualize the variation of the 
magnetic field components from the equations. In order to allow better 
visualization, the equations were programmed on an analog computer, for 
a satellite with an inertially fixed attitude, such as the 0A0. 

A portion of a typical record of the magnetic field in vehicle 
coordinates as a function of time for one orbit is shown in figure 1 for 
a near-earth circular orbit at an inclination of 35° to the equatorial 
plane. The vehicle coordinates are inertially fixed in space with the 
z axis parallel to the earth's polar axis. The magnitude of the magnetic 
field has been normalized with respect to B re;£ . which corresponds to the 
magnitude of the field intensity directly above the magnetic pole and is 
about 0. 25(l0 ) -4 webers/meter 2 at orbital altitude. Records of the type 
shown in figure 1 were taken for various orbital inclinations with the 
satellite rotating about the earth in the same direction the earth rotates 
about its own axis. The results in figure 2 show the variation over a 
24-hour period. As is shown in figure 1 the vector sum of the components 
varied slightly in magnitude while changing in direction. Equations (l) 
show that this variation cannot be more than 2:1 for a circular orbit. 

Each field component consists essentially of two sinusoids, one of which 
varies at approximately twice orbital frequency and the other at approxi- 
mately the frequency of earth rotation. These results can be visualized 
as follows: 
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Geogrophic axis 



Sketch (b) 


Sketch (a) shows the lines of force of the magnetic field, which 
are symmetrical about the magnetic pole. Equations (l) , along with the 
sketch, show that the magnitude of the field, at a given altitude, is 
a minimum on the magnetic equator. Sketch (b) shows that a satellite 
-will, in general, cross this equator twice per orbit. As a result, there 
will "be two minimums , and hence two cycles, per orbit. Because of the 
earth's rotation the period will only approximate twice orbital period. 

If the satellite orbits directly along the magnetic equator, the field 
at the satellite, in an inertially fixed vehicle coordinate system, will 
be constant. However, such an orbit can be achieved only momentarily 
since the magnetic axis moves in inertial space because of the earth's 
rotation. Therefore, for any orbit, the components of the magnetic field 
in an inerti all y fixed coordinate system can never remain constant. 
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As a polar orbit is approached, peculiar variations in the field 
are seen by the satellite. The magnetic field vector has its greatest 
magnitude at the magnetic poles. The resulting plot of the field 
components versus orbital time is a distorted sinusoidal wave as shown 
in figures 2(c) and (d) . 

It is apparent from this study that the prediction of the instanta- 
neous field is quite difficult, being a complicated function of time and 
the orbital elements. In addition, extraterrestrial effects, such as sun 
spots, will alter the magnetic field. Since knowledge of the magnetic- 
field vector is essential to the type of control system to be discussed, 
and since prediction of the field is unsatisfactory, the magnetic field 
must be measured aboard the satellite. 


TWO -AXIS CONTINUOUS CONTROL 


In the preceding section the earth’s magnetic field as a function 
of satellite position was discussed. In this section it will be .shown 
how this field can be used to generate a torque useful in controlling 
the attitude of the vehicle . Such a torque can be generated by develop- 
ing a magnetic moment, M, aboard the vehicle, which will interact with 
the earth's field. Accordingly, a known field will be postulated and a 
determination will be made of the magnetic moment which will give the 
desired values of torque. In the computation of the torque, the assump- 
tions will be made that there is a uniform, nonvarying field, B, and a 
constant moment, M, which have a constant relative orientation. For the 
nearly inertially fixed vehicle, the effects neglected, as a result of 
these assumptions, are small. 

Given a magnetic moment, M, and a magnetic_field, B, a torque will be 
created, equal to the vector or cross product, M X B (ref. 2). In order 
to control the magnitude and direction_of the torque one of these two 
quantities must be varied, and since M is developed aboard the satellite, 
it will be the quantity to be varied. For the purpose of discussion, we 
will assume that M is obtained in terms of its components, by passing 
currents through coils in the vehicle. Assume three coils, as shown in 
figure 3, which develop the three moments M x , My, M z along the vehicle 
axes, with 


M = iM x + jMy + kM z 
The torque developed is then 


T = M X B 


T = iTx + jTy + kT z 


where 
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and. 


B — iBx + jBy WBz 
The vector product can be expanded as 

T x := MyB z - M z By 

Ty = M Z B X _ M X B Z ^ 

T Z = M X By - MyB X 


( 2 ) 


The natural inclination is to equate the above torque components to 
the desired torque components, proportional to the angular errors about 
the three axe s . 


The three resulting equations, however, cannot be simultaneously 
satisfied for arbitrary values of desired torque because of the nature 
of the vector cross product. The resulting torque, by its definition, 
must lie in a plane perpendicular to B. This factor is one restriction 
on the vector, and only two additional restrictions may be specified. 
Thus, it is impossible to provide simultaneous proportional control about 
Hi 1 three axes by means of coils interacting with the magnetic field. 


For the particular case of interest here, namely, a vehicle with an 
inertially fixed attitude, whose x axis is to be pointed in a specified 
direction, Ty and T z should be directly commanded as a function of the 
y and z errors; that is, the desired torques are 


Ty = G(s)Ay 
T z = G(s)A z 


(3) 


where G(s) is that transfer function which will produce a satisfactory 
system response . Since the torque vector is then completely specified, 
equations (2) and (3) can be solved for the remaining torque component, 
T x , in terms of the field and the errors, giving 


T 


x 


-G(s) 


/^z®z + A y B y 
V . 


w 


This equation immediately points out a difficulty . If Ty and T z are 
always obtained as specified by equations (3) then at times when B x =0, 
T x = oo, which in turn means that T = oo and so M = which is not 
physically realizable. Thus, when B x = 0, or becomes sufficiently small, 
equations (3) can no longer be satisfied. For the moment, this condition 
will be neglected. 
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Of i S desired . to solve equations ( 2 ) and ( 3 ) for the three components 

re strict ion 6 on 1 ” ^M^can ^ here are ^ eally only two equations, one additional 
M cVin ij -k , an be made. One attractive restriction would he that 
uld he a minimum, that is, M should he perpendicular to B or 


M • B = 0 


(5) 


Another possible restriction, 
vehicle, is 


which only requires two coils aboard the 


A 

4 

7 

4 




The system developed from the first of these restrictions will he called 
the three-coil system, while that from the second will he called the two 

S&oT- Por the flr6t oase ' the - -3^ « 2 


If these 

to he 


Ty = M Z B X - M X B Z = G(s)A y 
T z = M x B y - M y B x = G(s)A z 
M X B X + MyBy + M Z B Z — 0 

equations are solved, the coil magnetic moments 


(7) 

are found 


M* = -G(s) 
My = -G(s) 
M z = G(s) 
for B x 0 


AyB z 

1b 2 


B 2 


A y B y B z + A z( B x 2 + B z 2 )' 


BvB 2 


B V B" 


A z ByB z Ay(B x 2 + By 2 ) 




B X B 2 


ByB 2 


C8) 


wLn h6 B 6 ^0 ^Len 11 ^ 3-+? US USedy Ve should consider what happens 

wnen B x 0 . When this condition occurs. My and M „ will become laree 

large" “o ter ' T? a ll5tlng SS aSf 

I“ ed iy Mx t raa °^?« ^ K ^ B,. these SrSs 5uTS±e 

system vhen* B ^ ^ d ) ^ J educe M x toward zero. Thus the contrS 
system, when B x - 0, reduces the quantity A V B 7 - A„B, r to zero Fnn 

e^mple, if in addition. By = 0, then Ay 7ill he oZtrol^Tto zero 
and A z will be uncontrolled. ’ 
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The second restriction, M x = 0 (two-coil system) , gives the following 
equations to he solved: 


T y » M Z B X = G(s)Ay 
T z = -MyB x = G(s)A z 


Solving these equations for M gives: 


M z = 




G(s)A y 

*T 

G( s ) A 2 


B- 


X 


for B x ^ 0 


( 9 ) 


( 10 ) 


These equations are simpler than the corresponding equations (8) for the 
three-coil system and require only knowledge of B x . In addition,, the 
controls for the two axes are uncoupled since each torque is a function 
of only one error . 

However , if these equations are used at times when B x = 0, both 
M z and My will be saturated and there will be no control of the vehicle 
about any axis. The two- and three-coil systems have similar limitations 
under the action of an externally applied torque . If such torque is 
perpendicular to B ; then the control system will exert a countertorque 
to hold the vehicle fixed; and the position error developed will be deter- 
mined by G(s) . If; however^ the external torque is applied parallel to 
B; the system will be unable to keep the vehicLe from rotating about an 
axis parallel to B. Since B is constantly changing in direction; the 
effect of the torque can be eliminated over a period of time . Further 
study is necessary to determine the short-term effects of this torque on 
the vehicle pointing errors . 

Thus both of these systems have difficulties when B x approaches 
zero and when an external torque is applied parallel to B. Further 
study is necessary to determine which of these systems; the three-coil 
system (using lower power; maintaining partial control when B x = 0; 
and with more complex computations); or the two-coil system (using more 
power; complete loss of control when B x = 0; but with simpler computation); 
is to be preferred for a particular application. The actual orbit of the 
vehicle must be considered in such a study; since this will control the 
times at which B x = 0. 
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THREE-AXIS INTERMITTENT CONTROL 


A 

k 

7 

k 


The previous section has shown that continuous control of the vehicle 
can he maintained magnetically about two axes, except for short periods 
when the field vector is perpendicular to the desired pointing axis . It 
is also possible to control all three axes intermittently but less 
accurately. However, for coarse control, or for the unloading of momentum 
stored in inertia wheels, the accuracy might well be good enough. 

Suppose it is desired to change the total momentum of the vehicle by 
an amount H, where H might be stored momentum to be dumped, or a desired 
momentum, as a function of error, such as to reduce the error to zero. 

If the equation for the magnetic moment of the coil is selected as 


_ K B X H) 

M = 

B 2 


(ID 


the torque developed will be 


T = M X B = - 


K(B X H) X B 
B 2 


= -KHj_ 


( 12 ) 


where H ± represents the component of H perpendicular to B. The 
resultant torque vector is parallel to Hj_ and thus cannot change that 
component of momentum parallel to the magnetic field, but will be such 
as to obtain the desired Hj_ . Since the magnetic field is continually 
changing its direction, the total momentum will be controlled. Thus 
current coils can provide a smooth control when periods of absence of 
torque about one axis over part of an orbit are of no consequence. (Gas 
jets can also be used for this purpose and do the job as well. However, 
the gas system necessitates carrying some form of compressed gas aboard 
and is thus limited in useful life , whereas the magnetic system can 
utilize solar energy and so can operate indefinitely.) This system will 
not be discussed further, although it appears to be a simple and straight- 
forward method for controlling, or dumping, vehicular angular momentum 
about all three axes. 


COIL CHARACTERISTICS 


II a magnetic moment is to be used, for torcjuing the vehicle ^ the 
coil weight and power required must be considered. The magnetic moment 
of a coil with no ferro-magnetic material of enclosed area a, current I 
and turns N is ’ 


M = NIa 


( 13 ) 
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where the mximuin magnitude of M required is determined by the particular 
satellite and its anticipated orbit. The maxirum magnitude of M is to 
be used in calculating the coil characteristics . Assume that the coil, 
has the following properties: a combined wire cross-section area, b, in 

meters 2 , wire length per turn, l , in meters, wire resistivity, p, in ohm- 
meters, and wire density, a, in kg/m 3 . Also let b be divided among N 
conductors each of cross-sectional area b/MN Then, 

R = (1^) 

b 


and, considering only the coil losses. 


P = I 2 R = 


(Nl) 2 Zp _ M 2 lp 

"h 


(15) 


If a and l are expressed as functions of the coil diameter 
(15) becomes 

_ 16M 2 

P ~ ltD 3 ! 


D, equation 


(16) 


In a similar manner mass can he expressed as, 

m ■= bed = jtDho 


( 17 ) 


It is considered desirable to keep both P and m as small as 
possible . Since their magnitudes are different functions of the same 
parameters, this can be accomplished if their product 

l6M 2 pc 


Pm = 


D 2 


(18) 


is kept small 


It can be seen that this product varies inversely as 

.—4 -1 . _ -t n 


and 


directly as M 2 , so that a large coil diameter and small maximum torque 
will tend to reduce the product. The coil area, b, may be used to control 
the trade between power and mass. The choice of b will thus depend, 
on the relative penalties associated with power and mass, although this 
choice may be restricted by considerations of coil heating. The available 
power supply voltage , E, will then be used to determine the number of 

turns , 


N = 


bED 

4Mp 


( 19 ) 


and may further restrict the choice of b. 
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Consider , as a numerical example , the case of the 0A0, where there 
will be disturbance torques due to gravity gradient, solar pressure, etc* 
An early estimate of this disturbance torque was 100 dyne cm (ref. 4). 
Later estimates have shown that the torques due to magnetic effects may 
be several times this value. In order for the system to recover from an 
initial error in a reasonable time, a torque of 5000 dyne cm is required. 
This value seems more than adequate to take care of the steady- state 
torques. If a magnetic field of 0.1 gauss is assumed, with a 2.44 meter 
diameter coil, only 10.7 ampere-turns are required to obtain the desired 
torque. With, for example, I = 100 ma, N = 107 turns of copper wire of 
gage Wo. 25 the power dissipated is 2.21 watts and the coil mass is 0.473 
kg. The power calculated for the particular case is quite small, and the 
A weight is reasonable . 

4 

7 

^ CONTROL SYSTEM RESPONSE 


The basic control system considered in both the analytical and 
experimental phases of the analysis is shown for a single axis in figure 
4 along with the over-all transfer function. The stability of the system 
is determined by the transfer functions as shown in the figure . The 
sensor provides a voltage proportional to the angular deviation from the 
desired line of sight to the star under observation. The signal voltage 
is amplified and added with its derivative to obtain the command torque . 

The magnetic moment is then varied by changing the coil currents (see 
the previously developed equation (8) or (10) ) so that the interaction 
of the earth* s magnetic field with the magnetic moment of the coil develops 
a torque which will reduce the error. Equations (7) and ( 9 ) both specify 
that the transfer function t/g(s)A is unity, even though the moments 
M are functions of B . 


The undamped natural frequency and dimensionless damping ratio, 
as determined from figure 4, are: 


L0 


n 




It can be seen that should be chosen to provide the desired degree 

of damping . 

A study of the magnetic control system that used three coils was 
made on an analog computer. For this study the constants of the system 
were selected to provide adequate damping and a reasonable settling time . 
Settling time is defined as the time from detection of the initial angle 
error to the time the error is reduced to within specified tracking 
accuracy. The prime concern of this analysis was to determine what 
effects changes in the earth 1 s magnetic field would have on the system 



response, and no attempt was made to find the optimum control system. 

The maximum current available for producing magnetic moment was assumed 
to be limited by an amplifier, so that the maximum torque available for 
a given magnetic field would be limited. To simulate this limitation, 
and also sensor saturation, limiters were incorporated in the computer 
program. Values of the field were determined from the data obtained 
from the field analysis previously described. (See fig. 2, for example.) 
During the initial phases of the study, the earth’s field was assumed 
constant during any one run and was changed only in direction between 
runs. This was acceptable since the field variation was small during 
any one run. System response to initial error and other conditions that 
the satellite may encounter, such as roll about the line of sight and 
disturbing steady-state torques acting on the satellite, were determined. 
Typical computer results are shown in figure 5> as plots of sensor output, 
in seconds of arc, against time. In figures 5(a), 5(b), and 6 the sensor 
output was limited to ±1.2 seconds of arc, while in figures 5(c) and (d) it 
was limited to ±3 seconds of arc. On these last two figures the ordinate 
scale is such that the initial portion of the transient is not shown. 

This change in limi t, value had very little effect on the general character- 
istics of the transient response. In each case the initial error was 
6 seconds of arc . 

Figure 5 ( a ) shows the system response when B x is fairly small. 

The response for larger values of B x is the sane. At this value of 
B x/B = 0 . 05 , however, the response of the system will be affected if 
either By or B z = 0, as shown in figure 5 (h) • Here By = 0 and the 
response of the z channel is somewhat worse, with more overshoot and 
a longer settling time . 

Figure 5(c) shows the effect of very much smaller B x . Here the 
response is worse, with a settling time of around 100 seconds. The 
response is still stable; however, figure 5(<i) snows the case where 
B x = 0, and the response is divergent. 

Computer tests were also conducted in which the rate of change of 
the field component was constant, which is a muci better approximation 
to the act ual conditions in orbit. This changing field did not affect 
the response of the system except when the satellite passed through a 
point at which B x was zero as shown in figure o. Instability will occur 
only when an angular error exists simultaneously with B x = 0. Due to 
the displacement of the earth's geomagnetic axis with respect to the 
polar axis, combined with the earth's rotation about the polar axis, the 
magnetic field components at the satellite in any practical orbit would 
always be changing . In particular , times at which B x is small enough 
to cause difficulties are few and of short duration; consequently, control 
can be maintained for a greater part of the orbit. 

Steady- state torques about the y or z axis will result in steady- 
state errors about the respective axis . These errors will be proportional 
to the torque necessary to cancel the existing disturbing torque. Also 
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such errors will couple torque about the line of sight. The magnitude 
of the roll torque depends upon the magnitude and direction of the 
angular errors, as can be seen from equation (T)» If angular rotation 
about the vehicle x axis is not to be tightly controlled, this torque 
may be tolerated. 

Errors can be made very small by increasing system gains or can be 
eliminated by inserting an integrating network. However, if an integrating 
network is used, a supplementary compensating network must be incorporated 
to maintain stability. 

An experimental model of a two-axis magnetic control system using 
two coils was built and tested on a platform supported by an air bearing. 

“t The results of the tests showed no significant deviation from the theoret- 

( ical analysis . 

l- 


SUMMARY OF RESULTS 


A study was made of the variation of the components of the earth's 
magnetic field at a satellite, in vehicle coordinates, for a vehicle with 
an inertially fixed orientation. These components vary in an appro xima tely 
sinusoidal fashion with time , but the exact variation is difficult to 
compute. Thus, on-board magnetometers must be used to determine the 
relative orientation of the magnetic field. 

Simultaneous proportional control about three axes cannot be achieved 
by the use of torque developed by the interaction of the earth's magnetic 
field with a vehicle-developed magnetic moment. Control is possible, 
however, about two axes with the third axis uncontrolled, except when the 
tracking line is perpendicular to the field, a condition which can only 
occur for short periods of time . 

The magnetic moment aboard the vehicle can be developed by means of 
current-carrying coils. If only two coils are used, the system is quite 
simple, but when the tracking line is perpendicular to the field, all 
control of the vehicle is lost . If three coils are used, the system is 
more complex, but less power is used and some control of the vehicle is 
retained when the field becomes perpendicular to the tracking line . In 
either case there will be an undesired torque about the tracking line, 
which may be acceptable . The power and weight requirements for the 
coils appear reasonable . 

It is possible to obtain torque about all three axes of the vehicle 
on an intermittent basis • This torque will change the vehicle angular 
momentum. The change in momentum can be used directly to produce a desired 
change of the vehicle orientation, or it may be used as a dump of stored 
momentum (desaturation of a wheel control system) . The equations for 
obtaining this type of control are relatively simple. 
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The use of the magnetic field for desaturation or fine control has 
the advantage that no fuel is required and the lifetime will be unlimited, 
discounting failures . In addition, for the fine control system, no auxil- 
iary system for momentum dumping is required. 

Ame6 Research Center 

Rational Aeronautics and Space Administration 
Moffett Field, Calif., May l8, 1961 
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appendix a 

DERIVATION OF THE EARTH’S MAGNETIC FIELD AS A FUNCTION OF 
SATELLITE COORDINATES 


it is I ?ir?tT lytiCal f tUdy ° f the effects of the earth’s magnetic field 
in sat! necessary to determine the components of the magnetic field ' 

o?bit hc ™ they ^ as £ L 

exuite in the satellite coordinate system. 

ss'H™ H; 
ss 

centered; end (4) a geographic coordimte^ySem^lich^e SrST^ntered. 

(fig.™) habits a V° U r ! “ ertial coordinate system 

Y v TI Zl P° inted at the north celestial vole anh 

X! - Yj define the equatorial plane. P ° le and 


maffneti 1 ^ ^S^tic coordinate system has its ^ axis along the earth’s 
magnetic dipole , positive toward the geographic north nnlp Qr ,n . 
oriented so that the satellite is almys ?n 7hJ Jt f 

with a negative component (fig. 8U pLane ’ 

g star, with Yg and Zg perpendicular (fig, 7). 

from the fLTSLr“Lt71 e te h (fig tB 9 ) XG TOCtOT 

y and . axes is interim in the follSg L‘“ YLSed. 

~ te 

rotations involved are shown in figure ^ 7 tvi^n^ n ^ ^ . . * j ‘* le 

about the appropriate axis, although in some cases a LgatirelnSris 0 " 16 
~ Is r fl^ Tfdu the bet “ e “ ihertia^and^vehicle 

eyate^TStl^rS/rr ^ ™ 
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System 2 Is then obtained by a further rotation about Yi Y 2 
through the angle -D s so that X 2 is alined with Xg. Thus RAg an 
Dc are the right ascension and declination of the point on the celestial 
sphere at which the satellite is looking. One additional rotation about 
X 2 = Xo through the angle cp s is necessary to completely define the 
vehicle coordinate system with respect to inertial space . 

Next, to define the magnetic system with respect to the inertial 
svstem. consider first the intermediate system 3 (fig- • I-kis 1S 
obtained by rotating about Zj = Z 3 through the angle RA M - 
RA M is the right ascension of the earth’s north magnetic pole and can 
be ^written as RA^ = R% p + w e^, where u> e is the angular velocity ^o 
earth rotation. System 4 is obtained from system 3 by a rotation about 
S =V 4 “tooSgh the angle -I* (equal to latitude of the earth' s north 
magnetic pole) so that X 4 is along %• A final rotation about 
X, = Xm to bring the satellite into the plane with a negative 

Z M is required . This angle is % and completes the specification of 
the orientation of the magnetic coordinate system with respect to the 
Inertial coordinate system. 

Finally, to define the geographic system with respect to the magnetic 
system, consider system 5 (fig- 8)- It is obtained by a r ° tatl ^ 

Ym = Y 5 through the angle until X 5 is aimed with Xg- The ang 
about X 5 = Xn between system 5 and the geographic coordinate system is 
not^r e quired • Thus, the angles RAm, L M , cpm, and ♦ define the location 
of the X axis of the geographic coordinate system, and hence locae 
the satellite with respect to inertial space. At this point % and t 
are -unknown. 

In matrix notation, if 1= I XG P is a vector P ointed to £ ard the 
satellite, its components in the inertial coorcinate system, Rj, are 
given by the following expression. 


A 

4 

7 

4 


Rj = [-RAmJ z [Lw] y [-9 m] x [-♦-‘y i-X&R 


(Al) 


The satellite can also be located with respect to inertial space through 
the orbital elements as follows (fig- 9 ) ' Sysuem 6 is obtained by a 
rotation about Zj = Z 6 through the angle «, alining the X 6 axis ng 
the line of nodes. Rotation about Xq = X 7 tlirough the inclinatio 
angle i locates system 7, with X 7 Y 7 defining the orbital plane . A 
further rotation through the angle 0 about Z 7 = Z 8 until X 8 is 
along Xg locates system 8. The angle 6 represents orbital motion, 


and can be written as 0-0 


;n as a - o D + w Q t * . 

The angle about X a = Xg between system 8 and the geographic coordinate 

system is not needed. 


where w 0 is orbital velocity 
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In matrix notation: 

R I = ^“ ± Jx T XG R (A 2) 

Equations (Al) and (A2) may now be equated and by premultiplying the 
resulting equation by [RA m ] and then [-L M ] the following expression 
is obtained: X 


I-%l x [-t] y 1 xg R = 1-Ly)^ [BA m ] z [-G] z [-i] x [-e] z i xg r (A3) 


A 

7 

4 


where the two sides_of the equation are different means for determining 
the components of R in coordinate system if. All angles on the right- 
hand side are known_and if the indicated multiplications are carried out, 
the components of R in coordinate system 4 can be determined. 


R = 


where R X4 , R y 4 > aQ l R Z4 represent 
system 4. 

• Explicitly carrying out the multiplication indicated on the left-hand 
side of equation (A3) gives 


R X4 

^4 

R Z 4 | 


(A4) 


the components of R in coordinate 


o 

o 

J 


c<Jr 0 sf 


1 ~ 



0 c<p -s9 


0 10 


0 

R = 

scpsi); 

_0 scp ccp 


_-si Jr 0 c vfr 


_ 0 _ 


_-CCpSl|/_ 


(A5) 


\ an< ^ S re P resen t cosine and sine, respectively. Equating (A4) 
and (A5J and solving for the cosine and sine of each of the angles gives 



scp 


/ 


% 4 

R v . 2 + R, 2 


-R* 




C(p = 


I 


R y/ + R z/ 


(A6) 


■V* ■ “Z 4 

so that. If desired, the angles cp and * can be explicitly determined. 
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Sketch (c) 


We have now located 
the satellite in the magnetic 
field by means of the angles 
cp and and the components 
of the magnetic field vector 
in satellite coordinates can 
now he determined. 

It is assumed that the 
magnetic field of the earth 
can he represented by a 
simple magnetic moment , with 
components as shown in the 
accompanying sketch. The 
axes X M and Z M of the mag- 
netic coordinate system are 
shown. The satellite is 
lccated by definition in the 
X #M Pl ane with a negative 
Zf.(. The field at the. satel- 
lite has the components 


2p 0 P H cos \|r 
Br = krtP 

p 0 P H sin * 

BT = W~ 


A 

4 

7 

4 


where Pg represents the magnetic moment. 

To determine the inertial components of B, the given polar components 
must first be rotated through the angle -t to the rectangular magnetic 
coordinates and then through the angles -<P, bn, and -RAM* Further 
rotations through the angles RA S , -D s , and <P S will then give the compo- 
nents of B in satellite coordinates resulting in the following equations. 


Bi = [-RA M ] z [L M ] [-cPm] x 


-2ci|f 

0 

s\|r 


^o P H 

4itR 3 


(A7) 


B s = [4> s ] x 


[-D s ] y [RA S ] Z Bi 


(a8) 


It should be noted that B in coordinate system 4 is a function of y 
and t only, and can be expressed in terms of the components of R m 
the same coordinate system- In particular; 
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B 4 = [-%] I-t] 

y 


-2c^ 

0 

st J 


4*R 3 


A. 

'+ 

7 

'+ 


-2c 2 t + s 2 ij/ 

■3scpctst 
3ccpctst 


^o p H 

4itR 3 


Substituting (A 6) into (A9) gives 


B, = - 


Then 


3R X . 


1 - 


R 2 
3Rx 4 R y 4 


R 2 

3 Rx 4 R Z 4 


R2 


M*: 


oH 


Ij-jtR 3 


B I = f-R%J z IlMly B 4 


(A9) 


(A10) 


(All) 


Combining equations (All), (AlO), and (A 3 ) then gives the component, 
of the earth's magnetic field, in satellite coordinates. 
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Figure 3.- Configuration of errors and coils in vehicle coor dina te system. 
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Open-loop response 
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TL = ■ 
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Closed loop response 
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Figure 4.- Block diagram of one axis of magnetic control system 
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Figure 5.- Continued. 
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(d) Bx/Bj-ef = 0, By/B^f = 0.36, and B z /Bj. ef = 0.51 
Figure 5 • “ Concluded • 




CvJ 


co <1- o 

i 


spuooas ojd 

‘" n *V 


00 CVJ CO ^ O 00 

r — ■ ii 

spuooas ojd 

« u,n *V 

ind*no josues 


CVJ 

I 


o 

6 



* 8 Jg 

*9 


— w 
o 

CVJ 


O 


i 


Figure 



36 



Figure 7*- Transformations between the inertial and the vehicle coordinate syste m s « 
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Figure 8.- Transformations between the inertial and geographic coordinate systems 

magnetic coordinate system. 








